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Abstract 



We have derived the radiative-transfer equation for a point source with a specified intensity and spec- 
trum, originating in the early Universe between the epochs of annihilation and recombination, at redshifts 
z B = 10 8 H-10 4 . The direct radiation of the source is separated from the diffuse radiation it produces. Optical 
distances from the source for Thomson scattering and bremsstrahlung absorption at the maximum of the 
thermal background radiation are calculated as a function of the redshift z.The distances grow sharply with 
decreasing z, approaching asymptotic values, the absorption distance increasing more slowly and reaching 
their limiting values at lower z. For the adopted z values, the optical parameters of the Universe can be 
described in a flat model with dusty material and radiation, and radiative transfer can be treated in a grey 
approximation. 

1 Formulation of the problem The action of the source causes the thermodynami- 



cal equilibrium in its vicinity to be violated, since the 
Let the early Universe be described by the standard source radiation is either nonequilibrium (for exam- 
model, i.e. by the Friedmann equations (see, for ex- pie, it has a power-law spectrum) or corresponds to 
ample, [1]). Let us suppose that a source of radiation a higher temperature. The radiation can experience 
is switched on at some time between the epochs of bremsstrahlung absorption and Thomson scattering, 
annihilation and recombination, that this source has In [2] , it was suggested that such sources could af- 
a specified geometry, spectrum, and total luminosity, feet the background radiation and hence be observed 
and that it radiates during some time interval. The via perturbations of this background, contrary to the 
intensity and spectrum of the source may be time- prevalent opinion that the energy of such sources is 
dependent. Such sources could occur in a number of totally dissipated. Both possibilities require careful 
processes, for example, during the amplification or verification using the methods of radiative-transfer 
damping of fluctuations in the matter distribution, theory. 

in the formation of primary black holes, etc. We assume that the source is pointlike. The 
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pointsourcc problem is also of independent interest: 
it is fundamental, since all other types of sources can 
be reduced to a set of point sources. 

If the source radiation is intense, it affects the sur- 
rounding matter in several ways. First, the radiation 
creates an additional pressure, which affects the lo- 
cal expansion of space. Second, this same pressure 
gives rise to inhomogeneity of the expansion. Third, 
a shock and, in the case of a periodic source, mat- 
ter oscillations may be formed. Fourth, the source 
can distort the metric of the surrounding space. In 
all these cases, the geometry of the source is the key 
factor. 

Here, we will assume that the source radiation is 
moderately intense and does not affect the expan- 
sion of the surrounding space. To start with, we will 
present the model for the Universe at the considered 
epoch. 



2 Two-component model of the 
Universe: dustlike matter 
and radiation 

2.1 BasicRelations 

We consider a model with a homogeneous and 
isotropic expanding Universe that contains equilib- 
rium plasma and Planck radiation in the stage when 
the plasma is radiation-dominated; more exactly, be- 
tween the epochs of annihilation (z < 10 8 ) and re- 
combination [z > 10 3 ). In this stage of the evolution 
of the Universe, the cosmological term does not play 
a marked role. 

At the temperatures corresponding to these red- 
shifts, the matter can be assumed to be nonrelativis- 
tic, and the conditions of thermodynamical balance 
with a common temperature for the matter and ra- 
diation arc fulfilled [1]. The pressure of the matter 
may be neglected, i.e., the matter can be considered 
to be dustlikc, in spite of its nearly total ionization. 

Under these conditions, in the usual notation, the 
evolution of the expansion is described by the equa- 
tion 
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-p u R 2 - kc 2 



(1) 



where a dot denotes a derivative with respect to the 
time t, the total mass density of the matter is equal to 
the sum of the matter and radiation densities, p u = 
pd + Pr, and the parameter k assumes the values 

1, for a closed model, 
0, for a flat model, 
— 1, for an open model. 

The matter density decreases in inverse proportion 
to the volume, while the expression for the radia- 
tion density contains another power of the radius of 
curvature in the denominator, resulting from the cos- 
mological redshift: 

n3 c>4 

n , - n a ill - EA n - - EjL 

p ~ Pd R 3 ~ " 3 ' Pr ~ Pr m ~ rA- 



R 4 



or 



Here and below, zero superscripts will denote cosmo- 
logical values at the current epoch. 

We will use the dimensionless time coordinate r] 
and dimensionless radius of curvature (scale factor) 
a(r]): 

R(t)dr, = cdt, o(rj) = ^. (2) 

Ho 

With these variables, Eq. (1) can be transformed into 
^{p> + P*)Rl-ka 2 . (3) 

This model was first considered by Chernin [3]. 

Let us also introduce the Hubble function and the 
critical density expressed in terms of this function: 

^ R a c 1 da c a' (77) 
R a a Rdi] Rq a 2 (77) ' 
_ 3ff 2 
Pc ~ 8^G' 

Along with the density itself, we will use the critical 
parameters 



fid = — , O r 



-, n u = n d + n r , (4) 

Pc Pc 

which can be used to rewrite (1) in another form: 

R 2 H 2 (l-n u ) =~kc 2 . (5) 
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2.2 Auxiliary Functions 



Let us consider the ratio 



Let us define several mutually related functions, in 
terms of which the solutions and relations of the 
model will be written. The two main functions are 
snfc(?7) and sck(r)), which are defined by the equalities 



sn fc (r;) 



sc fe (?7) 



sin rj when k = 1 , 

1] when k = 0, 

shr] when fc = — 1, 

1 — cos rj when k = 1 , 
r) 2 /2 when k = 0, 

ch 77 — 1 when fc = — 1 , 



(6) 



sc'fcO?) = sn fe (?7). 

Two functions are related to function (6). One of 
these, 

cos Tj when k = 1, 
cs fc (?7) = sn fc (?7) = { 1 when fc = 0, (7) 
ch?7 when fc = — 1, 

is its derivative, while the other, 

arcsiny, fc = 1, 

Wk{y) = { y, k = o, 

arshy = \n(y + yj\ + y 2 ), fc = -l, 



dar fc (y) 



dy y/l - ky 2 

is the inverse of snfe(x), i.e., arfc(snfc(?7)) = 
77, sn fe (ar fc (y))) = y. 
Another function, 

{rj — sini] when k = 1, 
rj 3 /6 when fc = 0, 
sh 77 — when fc = — 1, 

is the integral of scfc(77),i.e., cn' fe (??) = sck(rj). 

It can easily be verified that the following relations 
between the introduced functions are satisfied: 

cs2+fcsn| = l, (8) 
sn| +fcsc^ = 2 sc fe , cs fe +fcscfe = l. 



Tj 

ctg - when fc = 1, 

cth — when fc = —I. 
2 

Its derivative with respect to 77 assumes a simple form 
if we use the relation between the derivatives and the 
functions and equalities (8): 

sc fc -fcsc|-2sc fc +fcsc^ = _J_ 

SC^ SCfe 

Using the second of the equalities, this ratio can, in 
turn, be expressed in terms of y: 



dy y 2 + k 
d?] _ 2 ' 



(11) 



To obtain the correct solution, we must apply the 
appropriate initial condition: y = 00 when rj = O.The 
arguments of the functions in relations (8), (10), and 
(11) have been omitted for brevity. 

2.3 Solution of the Equation of Mo- 
tion 

It was shown in [3] that the solution of (3) in implicit 
form is 



a{rj) = a r sn k (rj) + a d sc k {v), 
R 

t= — [a r sc fe (rj) + a d cn k (rj)] . 



(12) 



Indeed, we can easily verify using relations (8) that 
the function a(rj) satisfies (3) if we assume 
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The obtained coefficients are expressed in terms of 
the critical parameters. Let us first derive an expres- 
sion for the current radius of curvature [1]. To this 
end, we apply formula (5) to the current epoch: 



rt a K 



kc 2 



HqRq 
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We now have, in accordance with the definitions (4), 
4ttG 



ad 



3c 2 



(i o R 2 _ 47TG 3H§ 2 



k n° d 



3c 2 8ttG 



3c 2 



21-ns 2|i -nsr 
p"q,^Rq = - n - — ^n^Rn 



,2_ 8^^ 0n0D2 _ 8^G3H 2 

3c 2 87rG W " 



i-fiS |i-fig|' 

The inverse formulas have the form: 



= 2a d 

d 2a d + a 2 -fc' 



0' = 



2ad + a 2 — k ' 



0° = 



2a d + a 2 
2ad + a 2 — k 



In the following subsections (to Subsection 9), we will 
assume that k ^ 0, although some functions are cal- 
culated using formulas that include the special case 
of a flat space. 

2.4 Relation to the Redshift 

Let us derive expressions for the quantities associated 
with the solution (12) in terms of the redshift z. By 

definition a(n) = — — = . Successively substi- 

Ro 1 + z 

tuting in (12) the expression for sn^ or sc^ in terms 
of the other, and also using the last relation in (8), 
we obtain 



Q-B(l-kA) 
1 + kB 2 ' 
B 2 + A - BQ 
1 + kB 2 ' 
l-kA + kBQ 
1 + kB 2 ' 

Here, we have introduced the notation 
1 1 2|1-Q°| 1 



sn fc (ry) = 
sc k (ij) = 
cs k (v) = 



a d 1 + z 



1 + z' 



B= a L = 2^/n°\i-n°\ 



ad 



(13) 



Q = x/B 2 + 2A-kA 2 . 



The derivative a'(rj) and the Hubble function can be 
expressed in terms of the redshift in this same way: 

a'{rj) = a r cs fc (?7) + a d sc k {ri), 
H c a'(r]) a'{r]) 



H H R a 2 ( V ) a 2 ( V ) 



The ratio (9) and the product of this ratio with B 
can also be expressed in terms of z: 



y(v) = 

X = By = 



kAB -B + Q 
B 2 + A - BQ ' 
kAB -B + Q 
B + A/B-Q- 



The critical parameters at the current epoch (B is 
constant) can likewise be expressed in terms of the 
quantities in (14): 



2A Q 

QV 

2A + B 2 



Ql 



A 2 

1 - Q° - — k — 

It is easy to verify that the identities a(rj ) — 1 and 
a '(Vo) = HqRq/c = l/\/|l — ^21 are satisfied. 

2.5 SpecificModels 

According to the above model, the curvature of space 
can be arbitrary (constant). However, recent obser- 
vations lead to the conclusion that space was flat over 
the entire history of the Universe. In this context, we 
will consider two specific models. One is open, close 
to flat, and does not take into account the vacuum 
component (we will call this the nonvacuum model), 
while the other has three components and is strictly 
flat, but with the vacuum neglected at z > 10 3 . 

We will adopt for the Hubble constant H = 
65 kms -1 Mpc -1 , the corresponding critical density 
p° c = 7.940 x 1CT 30 g/cm 3 , and the critical parame- 
ters for the baryon component = 0.025 and dust- 
like component fl d = 0.25. The baryon density is 
then pi = 4.684 x 10~ 31 g/cm 3 , the density of dust- 
like matter is p° d = 4.684 x 10~ 30 g/cm 3 , the mass 
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density of the radiation corresponding to the temper- 
ature T = 2.7277 K, is p° = 4.63 x 1(T 34 g/cm 3 , and 
fi? = 5.83 x l(r 5 . 



2.5.1 Nonvacuum model. 



-n° = i-fi§- 



For the adopted model parameters, 1 

= 0.7499; i.e., this model is open (k = — 1), which 
is natural, since there is no vacuum component. The 
other parameters are a r — 8.869 x and ad = 

0.1667. 

2.5.2 Strictly flat model. 

For k = ,we also have 1 — 51 u = 0, so that many of 
the introduced values become equal to zero or infinity. 
In this case, the limit transition is complicated, and 
it is easier to consider this case individually. 

We will assume that the most plausible model of 
the Universe for redshifts less than 10-100, in partic- 
ular for the current epoch, is a flat model with three 
noninteracting components: dustlike matter, radia- 
tion, and the vacuum. Based on this model, we will 
also construct a two-component flat model for the 
considered epoch of radiation-dominated plasma in 
which the two components have the same densities, 
but there is no vacuum component. All values cor- 
responding to this model will be denoted by a tilda. 
The contribution of the vacuum can be neglected for 
the redshifts considered (10 3 < z < 10 8 ), so that the 
Hubble functions in the two-component and three- 
component flat models differ only slightly: 



8ttG 



u = w^v^i + ^ + ^i + ^ + pa, 

As 

H -- 



8irG 



p°(l + z) 3 + p?(l + z) 4 . 



We will use the same critical parameters and 
as above; the critical parameter for the vacuum is 
then fi A = 1 — — Q° ■ In the adopted model, these 
values are: 



H n — Hn 



P? 



P° d + P?+P°A 



iWi 



(14) 
(15) 



ng 



i 



n? 



i 



"A 



(16) 



For brevity, let us set Cl^ = CI, in which case f2 r = 
1 - O, so that tt^ = Q and 0° = 1 - Q . 

The radius of curvature R and its current value i?o 
are meaningless in a flat Universe (they are infinite, 
and can be selected arbitrarily in the formulas), and 
should not arise in expressions for physical values. 
Indeed, the solution (12) can be written in the form 



a(rj) = a r ?? + a d 



-—HoRo 
\l rj 



where we have introduced a new time variable that 
is linearly related to the previous time variable and 
is expressed in terms of a(ri) and the redshift: 

HqRq a 

C=^ — v 



2c 



3.1 



y/i-n (i-a) + VT=T%; 
i 

Radiative-transfer equation 
in the early Universe 

General Form of the Radiative- 
Transfer Equation for a Point 
Source 



Let us assume the source to be pointlike and 
isotropic. In this case, the problem is spherically sym- 
metrical, so that, generally speaking, the approach of 
the Tolman-Bondi model should be used. However, 
to first approximation, we can assume that the source 
radiation does not affect the metric of space and its 
expansion, and the evolution of the source radiation 
can be considered against the background of the stan- 
dard cosmological two-component model described in 
Section 2. 

We will describe the radiation with the mean oc- 
cupation number of photon states n, rather then 
the intensity. The advantage of this quantity is 
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that it is dimensionlcss and relativistically invariant 
(scalar). Due to the symmetry of our problem, the 
mean occupation number of photon states will de- 
pend on the time t or time variable 77, the distance of 
the point considered from the coordinate origin (de- 
scribed by the parameter \), the angle 9 between the 
radial direction and the ray along which the radia- 
tion propagates, and the dimcnsionless frequency x: 
n = n(rj, x, 6, x). The radiative-transfer equation has 
the form 

dn dn dn dx dnd\ dn dO 

di = ~dt + d^dt + d^~dt + 'd9 at " c c ' 

where I c is the collision integral, i.e., the difference 
between the numbers of photons with specified pa- 
rameters entering and leaving the state. 

We will use the parameter r\ instead of the time; 
the derivatives will be calculated with respect to r\ 
using relation (2), so that 



d_ 

di 



d 



r( v ) dr]' 



The derivative of the frequency with respect to the 
time is the easiest to find. Due to the redshift, 



v = vo 



Rq _ V _ Vq Rq 

R ' v± R 

t = - XoRo i l z = - xH > 



x - 



Ro 
~R' 



R 

where H = — is the Hubble function and v+ is some 
R 

specified frequency. 

To determine the other derivatives, we must ana- 
lyze the variations of the corresponding values along 
the ray. 

3.1.1 Left-hand Side of the Transfer Equa- 
tion 

Let us consider the ray along which the radiation 
propagates, having, in general, already been scat- 
tered. An arbitrary ray can be described using sev- 
eral parameters, one of which can be chosen to be the 
coordinates of the point on the ray closest to a point 
source at the coordinate origin. We will call this the 



initial point. Due to the symmetry of the problem, it 
is sufficient to specify this distance using the coordi- 
nate Xo- Let a photon pass through the initial point 
at time t Q — t(r) ). Its position on the ray will be 
described by the parameter x* associated with this 
time and measured from time t . Let us determine 
the derivatives with respect to x* of the variables 
describing the position of the photon on the ray. 




Fig. 1: The curvilinear triangle. 

To this end, let us consider the curvilinear triangle 
presented in Fig. 1, with its sides extending from the 
source S to the initial point, from the source to the 
photon at time 77, and from the initial point to the 
photon. Since all the lengths are proportional to the 
radius of curvature, we will consider dimcnsionless 
distances described by the coordinates xo, an d X*-> 
respectively; the corresponding dimcnsionless lengths 
are a = sn fe (x Q ), sn fe (x), and sn k (x*)- In the same 
order, the angles opposing these sides are those be- 
tween the ray and radial direction, 8, between the 
radial direction toward the initial point and the ray 
(a right angle), and between the radial direction to- 
ward the initial point and the direction toward the 
photon's position, 9*. 

We can write the geometrical integrals of the trans- 
fer equation for the constructed right-angle curvilin- 
ear triangle: 



= sn fe (xo) = sn fe (x) sin#, 
snfe(xo) 



CSfc(Xo) 



CSfe(Xo) 

= sn fc (x) 
sn fe (x*) 



= sn fc (x*)tg 6, 



cosO- 



sm( 



(18) 
(19) 

(20) 
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csfe(x) = r cs fe (x*), 

sn fc(x) a sn fe(x*) 
■ cos 6 = 



cos( 



csfe(x*) ' 
cs fc (x*)sin6». 



The functions sn/t and csfc are specified by equalities 
(6) and (7). Relation (20) follows from the two pre- 
ceding expressions (18) and (19). These relations are 
trivial in the case of a flat space. For a space with 
positive curvature, the triangle can be placed on a 
unit sphere in three-dimensional space, with the sub- 
sequent use of spherical geometry. To make the trans- 
lation to the case of negative curvature, the trigono- 
metric functions expressing the sides of the triangle 
must be changed to hyperbolic functions. 

The following derivatives can easily be obtained 
from these relations: 



4^ = cos e, 
dx* 

d0 ■ fl cs fc (x) 
- — = - sin 6» -—, 

dx* sn fc (x) 

dfl* sin 9 



(21) 

(22) 

(23) 
dx* sn fe (x)' 

The equation of motion of a photon along a ray is 
X* = t) — rj , so that the derivatives (21), (22), and 
(23) are simultaneously the derivatives with respect 
to the time coordinate r\. 

Thus, the transfer equation (17) takes the form 



dn dn 
— + cos9— 
dn dx 

- X H dn 



? cs fc (x) dn 
sn k (x) d9 



tile. 



dx 

For a flat space, when k = 0, sno(x) = X-, 
cs o(x) = 1) tne equation takes the usual form. 



(24) 
and 



3.2 Right-hand Side of the Transfer 
Equation 

In the problem we are considering, bremsstrahlung 
absorption and emission, as well as Thomson scat- 
tering, are important processes in the interaction be- 
tween the matter and point-source radiation. 



Bremsstrahlung absorption by thermal nonrcla- 
tivistic electrons with a Maxwellian velocity distri- 
bution and the temperature T is described by the 
formula [4] 

"cc = n e n + k cc (u, T), 
a cc = (l-e- h "/ k * T )a* cc , 



where the cross section is 



ky.,. 



fc" =8tt 



cTi (67rmfcB) 3 / 2 



k cc {v,T) - T1/ C 2 C ^ 3 , 
= 3.69xl0 8 cm 5 K 1 /^ 3 . 



The absorption coefficient has been corrected for 
stimulated emission. The Gaunt factor is taken to be 
unity. The bremsstrahlung coefficient in the occupa- 
tion numbers is related to the absorption coefficient 
by the Kirchhoff-Planck relation: 



^cc — Q^cc^ 



-hv/k B T 



(25) 



The absorption coefficient associated with Thomson 
scattering is ax = n c a^, where the cross section is 
e 2 



8^ / ~ 2 ^ 2 



3 V mc 



6.66 x 10" 



The indi- 



catrix for Thomson scattering is the Raylcigh index. 
Since the radiation field does not depend on the az- 
imuth, the equation will contain the azimuthaveraged 
indicatrix 

p{6,0') = l + ^P 2 (cos0)P 2 (cos0'), 

where i-2(/u) is the second Legendre polynomial. 

Thus, the right-hand side of the transfer equation 
has the form 

J c = -[a cc (u,T) + a T ]n(r],x,0,x) 
+e cc (v,T) + a T I + e*(n, X ,6,x), (26) 



n 

1 J sm6'p(6,6')&6'n(n,x,0', 



where describes the primary radiation from the 
point source. 
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The integral in (26) can be divided into two terms frequencies) and e* is 
containing two angular moments n(r), \, 6, x): 



no 



(n,x,x) = ^ J sm6'd6'n(r),x,0',x), 
o 



L = 



J dt J d 3 r J x 3 dxe*(r),x,0,x). 



The factor before the integrals is introduced to trans- 
late the occupation number into intensity. Express- 
ing time in terms of ij and writing the integral 
u over space in terms of integrals over the three vari- 

Let us now specify the source function e* and sep- ab , les ' ta } [n S into account the volume clement d 3 r = 



«2(»?, X, x) =- lsm6'P 2 (cos e')d6'n{r}, x, 0', x). 



arate the direct from the diffuse radiation. 



3.3 Direct Radiation 

The entire space is filled by homogeneous and 
isotropic equilibrium cosmological radiation with the 
mean occupation number 



n e {v,x) = 



e /iL/,x/(fe B T) _ ^ hi/*x /k B T _ j ' 



(27) 



It is straightforward to verify that this function satis- 
fies (24) with the right-hand side (26) without the last 
term. In fact, both the left-hand and right-hand sides 
vanish, the left-hand side because the mean occupa- 
tion number for the background radiation does not 
depend on anything, and the right-hand side due to 
the Kirchhoff-Planck relation (25). This is the equi- 
librium background radiation. In spite of the fact 
that it interacts with matter, overall, the matter and 
radiation do not affect one another during this inter- 
action; only a common temperature is established. 
After the epoch of recombination, the electron den- 
sity becomes extremely small (it can be assumed to 
be zero), and the right-hand side of the transfer equa- 
tion no longer contains terms describing the interac- 
tion between radiation and matter (the radiation is 
separated from the matter), while the left-hand side 
corresponds to the free propagation of the radiation. 
The quantity (27) remains constant and satisfies the 
transfer equation. 

Let us now consider the direct radiation of the 
source. We will first find an expression for e*. The 
relation between the total energy L emitted by the 
source (for all time, in all space, in all directions and 



R 3 (r])sn1(x)sm9dxd9dip, we obtain 



oo oo 

L = 2tt^ J R 4 (v)dv J snt(x)d X 



7T OO 

x J sintfdfl J x 3 dxe*(r],x,0,x). (28) 



It follows from the expression for the right-hand side 
of the transfer equation (26) that the units of e* are 
inverse length. 

If the power of the source radiation at the time t s 
with the time coordinate rj s at frequency x s is L6(t — 
i s ), we can write 



e*(ri,x,0,x) = L 



5{n - ?7 S ) 5(x) 



2hvi i? 4 (r7s) sn 2 fe ( x ) 
(5(cos# — 1) S(x — x s ) 
2ir xf 

and relation (28) is satisfied. 

Further, let us determine the mean occupation 
number for the direct radiation from the source. It is 
specified by the transfer equation (24) with the right- 
hand side (26) without any terms describing scatter- 
ing: 



n d 

— + cos6>— 
on ox 



sin i 



, cs fc (x) d 
sn fe ( X ) 80 



a(rj) dx 

= -[n + k cc (v,T) + a T ]n c R(iq)n, : {-q,x-,0-,x) 
+e*{r],x,6,x)R(r)). 
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Direct substitution verifies that the solution of this 
equation is given by the function 

■ A 6(x-V + Vs) @(x) 



n*(r],x,0,x) = L 
S(xa(n)/a(r] s 



2hut R 3 (Vs) sn 2 fe ( X ) 
) - x s ) (5(cos6> - 1) _ T 



x^ 2ir 

where the optical distance from the source to the 
point is (see below) 

t = t(t], t] - i] s , x s a(r] s )/a(r))) = r T (r?) - T T (rj s ) 

k%[E( V ) - E(rk)] 



3.4 



^xlT s 1/2 aV^ Vs ) 



Discriminating between Direct 
and Diffuse Radiation 



Let us represent the total occupation number for the 
photon states created by the source as the sum of two 
terms: 

n(v,X,0,x) = n*(r),x,0,x) + n d (n, \, 9, x), 

where rid describes the diffuse part of the radiation. 
This function satisfies the same transfer equation (24) 
as the initial function, but with the "source" term e* 
replaced by e s (''7, X, 9, x) — the part of the emission 
coefficient due to scattering of the direct radiation 
from the point source. This quantity is specified by 
the formula 



e s (r),x,0,x) = a T 



l + ip 2 (cos0) 



where the zeroth moment of the direct intensity is 



L c 3 Six-y + r/s) 9(x) 



^S(xa(T])/a(f] s ) -x s ) ^_ 

X r G 



The second moment of the direct radiation exceeds 
the zeroth moment by a factor of five, since P2 (1) = 1. 

For a noninstantancous and nonmonochromatic 
source, the emitted energy may depend on % and 
x s , and all the equations and their solutions must 
be integrated over these parameters. This does not 
present a problem, due to the presence of S functions 
in these expressions. 



3.5 A Planar Source 

For comparison, a layered planar source can be con- 
sidered. This is only possible if the geometry of the 
space itself is flat, fc = 0. As always in a planar ge- 
ometry, the equation of radiative transfer does not 
contain an angular derivative, since the angle 9 be- 
tween the direction of propagation of the radiation 
and the normal to the source plane along the ray 
remains constant. Thus, the transfer equations for 
a planar geometry can be obtained from the corre- 
sponding equations for spherical geometry, simply by 
omitting the term with the derivative with respect to 
the angle 9. In the expression for the source term and 
the intensity of the direct radiation, the combination 
<5(x)/sn|(ry) must be replaced by 5(x — Xs)- 
Thus, the transfer equation acquires the form 



dn n dn 
— +costf— 
or] dx. 



dn 

xH lT 
ox 



cX, 



where the collision integral l c is specified by the same 
formula (26). In this case, the source term is 



e*(r],x,IJ;x) = L. 



c 3 5(t]-t] s ) 



Hx-x s ) 



2hvt R 2 {n. 
S(cosr] — cos?9 s ) S(x — x s ) 
2ir x 3 

The source is assumed to be located at the level with 
Xs- Here, we calculate its luminosity per area of the 
planar boundary; for this reason, the denominator 
contains a second power of the radius of curvature. 
We assume that the source radiation does not depend 
on azimuth and propagates at an angle i) s relative to 
a normal to the layers. 



4 Optical distances 

4.1 Optical Distance from the Source 

Under the adopted conditions, the densities of elec- 
trons n c and protons n + depend only on time (and, 
of course, the cosmological model). In the case of 
total ionization, n + = n c = n°/a 3 (77); in the case 
of partial ionization, the ionization equation must be 
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solved. We will assume that the matter consists of 
completely ionized hydrogen. 

Expressions for the coefficient of bremsstrahlung 
absorption by thermal nonrclativistic electrons and 
the absorption coefficient for Thomson scattering 
were given in Section 2. 

Since the scattering coefficient docs not depend 
on frequency, the corresponding optical distance can 
be determined very easily. If a photon was emitted 
from the source at time i] s , its equation of motion is 
X = V ~ Vs- In the course of its motion, the opti- 
cal distance between the photon and the source will 
increase according to the relation 

x 

tt(v,Vs) = J aT(Vs+x )R(Vs + x )dx 



where the universal Thomson opacity is 



coefficient is a known function of 77 or z: 



v 

t t(v) = J a T (rj)R(rj)dr] 

Vi 

■n 

= <7 T J n c (rj)R(rj)dr]. 



(29) 



Since only differences in the optical distances are of 
interest, the lower integration limit in this formula 
can be chosen arbitrarily. 

The situation with the bremsstrahlung mechanism 
is more complicated, since in this case, the absorp- 
tion cross section depends on both the frequency 
and temperature. The time dependence of the fre- 
quency is described as follows. Let the frequency of 
a photon emitted from the source at time t s = t(rj s ) 
be v s . Then, its frequency at time t = t{r}), will 
be v — v s a(r) s ) / 'a(rj) . The temperature also de- 
pends on the time, according to the usual formula, 
T = T /a{r])=T s a(r] s )/a(rj), where T s = T M%) is 
the temperature when the photon is emitted. The 
ratio v/T = v s /T s is time independent. 

Since the particle densities, frequency, and temper- 
ature depend on time in a known way, the absorption 



: (»?,?7 B ,a;)=( 



1-e 



-hv,x s /knT t 



X 

) jn c {rj s + X ') 



^ + (VS + X') 3 3 7 . Tl/2 « V2 (VS + X') 

vfx^a 7 / 2 (r] s )T s ' 
xR(r, s + X ')dx'= (l - e-^'AB^ Rq 
xk cc {iy ir x s a(%),T s a(r] s ))[E(r]) - E(rj s )}, 



E(rf) = I n c (?y)n + (r/)a 9/2 (?7)d?7 (30) 



where 



is a quantity similar to the emission measure. The 
total optical distance from the source is 

r(r],r] s ,x) = r T (??,r? s ) + r cc (r], rj s , x). 

4.2 Calculation of the Scattering Op- 
tical Distance 

Let us choose for the initial times in (29) and (30) the 

times t; = t{rj\) for which the ratio y = snk(r])/ scfc(?y) 

assumes the values from the table. The table also 

a r 

presents the products X = — y = By for the same 

ad . 

times. Here, we have chosen simple values for this 
ratio; in fact, the values of r\\ are maximum for k = 
and fc = — 1, 7/ < so that the optical distance (29) 
will be negative. This does not cause any problems, 
since only the differences of these quantities, which 
are positive, will be needed. When k = 1, the value 
of rji is finite, since the range of variations of this 
variable is finite in a closed model. 

We will rewrite the integral (29), substituting 
the time dependence of the electron density, n c = 
n e/° 3 ( 7 ?) = (Pb/ TO H)/a 3 (??), where p° is the current 
baryon density: 





T T (rf) =a T n c R(T])dT] = <j t -^R / - 
J ma J a 



drj 



'(V)' 



Introducing the new variable of integration X = 

—y = By, where y is specified by (9), and using 
CM 
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Table 1: f]i,yi and Xi 



k 


v 


m 


Vi 


X; 


k = 1 


ctg(r?/2) 


7T 








k = 


2/77 


OO 








k = -1 


cth(r?/2) 


oo 


1 


a r /a d 



the characteristics of this ratio enables us to reduce 
the integral to the form 

/ dr? _1 f y 2 + k 

J a 2 (77) 2 7 (a r2 / + a d ) 2 V 
m yi 

= ~^MX,B,k), 

Calculating the integral within the limits indicated 
in the table yields for k = and 1: 

X 4- 2 + kR 2 

f T (X,B,k)=X—-^- 21n(X + l). (31) 

A.+L 

When k = -1, the constant B(B - 2) + 2 ln(l + B) is 
added, so that, in an open model, this function can 
also be written in a form that explicitly vanishes for 
the limiting value of the argument, X = B: 

h(x,B, i: :.v b) X ' u 



X + l 



-2 In 



X + l 
B + l' 



0, 



For arguments close to the limiting value, X 
the function (31) is equivalent to Jt{X, B, k) <~ 
kB 2 X(l - X) + (1/3 + kB 2 )X 3 , while for an open 
model, we obtain as X — > B then fi(X, B, — 1) <~ 



{x-by 
{B + iy 



B + 



1/3 -B 



(X-B) 



B + l 

Thus, the optical distance based on Thomson scat- 
tering is given by the formula 



a T p° h p° d /SttG „ 



4 m H pj) V 



Pr 



-1/2 



pj? \ 3c 2 
x[/ T (X s ,B,fc)-/ T (X,B,fc)] 



a,. 



where B = — , X = By 



B^44, and X < X s 



a d sc fe (?7) : 
for f] > r] s . The factor before the bracket can also be 

written as ^A^JL 



4 m H (ft°) 3/2 H ' 

4.3 Calculation of the Absorption 
Optical Distance 

Let us calculate the function (30). After making the 
same substitutions as for the Thomson-scattering dis- 
tance calculation, we obtain 



A 

m H 



A) f d7 i 

run) J a 3 / 2 {r/) 

Vi 

2 y i 

' i r ^RFTk 



V2j (a r y + a d ) 3 / 2 

Vi 



dy 



A 



1 a 



1/2 



V2 a? 



-f c (X,B,k), 



where 



x 



.UX,B,k) 



y/X 2 + kB 2 

{X l + 1) 3 / 2 



dXi. 



Only for B = 0(i.e. for k = 0) can the integral be 
expressed in terms of elementary functions: 

/ C (X) = / C (X,0,0) = 2VA 7 TT 



Vx+i 

2 



-4 = 2 



(Vx+T-i) 2 
x 2 



Vx+i(vx+i+iy 



(33) 



If B 7^ 0, the integral can be expressed in terms of 
elliptical integrals, but it is easier to calculate it nu- 
merically. 

4.4 Calculations of Optical Distances 
for the Nonvacuum Model 



(32) 



For the values in (32), we will obtain B 
1(T 2 . 



5.321 x 
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Let us consider optical distances in the selected 
model. We will first calculate the dimensionless func- 
tions fx and f c for redshifts z from 10 s to 10 2 . As- 
suming that the photon was emitted from the source 
at time rj s , the Thomson-scattering distance at time 
77 is specified by (32). The factor before the bracket 
in this formula is equal to 156.2. 

The optical distance based on absorption can be 
determined in a similar way. After substituting ex- 
pressions for its terms, the product RoE{rj) is shown 
to be 



but, since B = and k = 0, it does not contain any 
parameters: 

h(X)=X^±-2\a{X + \). 

Taking into account the relations between the Hubble 
constant and critical parameters (14), (15), and (16), 
we will see that the coefficients before the brackets in 
(35) and (32) also exactly coincide. 

The absorption distance is calculated in exactly the 
same way. This product is equal to 



R E{ V ) = 



1 



Pb 
m H 



8nG 



3c 2 



-1/2 



R E(r]) 



1 



xf c (X,B,-l) = -U4-) 



W d c 



P°b 



'fin 



2 \ m H / 1 - fin H 



~rfc(X), 



2 \rnnj fio H 
xf c (X,B,-l). 

Let us estimate the optical distance for the frequency 
at which the Planck function reaches its maximum 
at a certain cosmological time, v ra = c\^k^T /h, cw = 
2.821438. For this distance, 



where the function f c (X) is specified by the previous 
formula (33). The resulting formula for the absorp- 
tion distance is 



= (1 - e~ Cw ) 



(A 



^[ fc (x s ) - f c (x)}. 



(36) 



: (1 



c w \ I Pb 



TV 2 \cwk B J fi° 



1 - fio 2H 

Here, the coefficients in (34) and (36) also coincide. 



2H< 



[f c {X s ,B,-l)-f c (X,B,-l)}. (34) 4 6 Comparing the Two Models 



4.5 



Calculation of Optical Distance 
for a Flat Model 



The optical distances needed for our model have al- 
ready been determined. Let us rewrite them using 
new notation. The scattering optical distance is 



The optical distances in the two models considered 
differ only in the functions /;in the flat model, B = 0, 
while this parameter is B « 0.05 in the open model, 
so that the difference is insignificant. Indeed, this 
difference 



or no— 

to h H a R 



f T (X)-f T (X,B,-l) 



B 



X 



fin 



X + l 



yi -fi 2 

fio Cs 

_ ctt pI 



h 



(1 - fio) 3 / 2 
/ A /T~fi ' 2 
I fio C 
fio 



4 m H Ho (1 - fio) 3/2 

x[f T (x s ) - f T (x)]. 



(35) 



We have obtained the same formula, (32); the func- 
tion fx coincides with that determined from (31), 



-B(B- 2) - 21n(B + 1) 

is small for both small and large X. The term that 
is independent of X disappears when the difference 
of the functions is taken. The difference between the 
functions f c is more difficult to estimate, since these 
functions have different domains: 

B 

f c (X)-f c (X,B,-l) = J — ^373 dXa 
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x 



+ 



J {Xi + 1; 



X? - B 2 



)3/2 
B 2 



-dXi 



+B " l (X 1 + 1)3/2 Xi 



The order of the difference is £? 2 . For small X <~ 
B, this can be seen directly, while, in the case of 
larger X, the contribution of large values of X\ to 
the integral is small. 

We can describe the behavior of our functions for 
small and large X. When X — > 0, 



4-1 

/ T (X)=X 3 ^(-1)"^X", 

n=0 



fc {x) = x*f{-ir {2 l +1) \ . 

J y ' ^ y ' 2"n!(n + 2) 



n=0 



For large X, 



X ^ X™ n + 1 : 

n=0 

f c (X) = 2X 1 / 2 V- (_i)»-i ( 2 "-3)' 2n+1 _ 4 . 
J v ; ^ V y (2n)! X" 

n— v ' 

Here, as usual, we assume that (—3)! = —1, (—1)! = 
1. 

For the large redshifts of interest for us or, more 
exactly, for (1 — Qo)z ^ 1> the parameters £, y, and 
X behave as follows: 



c 



1 1 



y = ^ ~ 2 VI - £V, 

VT~T^ , 1-Oo 
A = 2 y ~ 4 — z. 



n 



fl 



Accordingly, the functions / are easy to determine 
for such z values. The scattering optical distances 
calculated for the open and flat models are only dif- 
ferent for the largest z, and even then only slightly: 



for z = 10 8 , in the fourth significant figure. The ab- 
sorption distances display larger but also insignificant 
differences, in the third significant figure for all z. 

Figure 2 presents the results of these calculations. 
The starting value for the curves is z = z s ; lgz s = 
n s = 4(1)8. To present the curves on comparable 
scales, the tt values on the curve with lg z s = n s 
have been multiplied by 10 1_n % while the r cc curves 
have been multiplied by 10( 17_ ™ s )/ 2 .The difference 
between these factors reflects the fact that the ab- 
sorption coefficient is proportional to T -1 / 2 . 




Fig. 2: Scattering and absorption optical distances 
as a function of the photon rcdshift in the course of 
the photon's motion away from the source. The tt 
values are multiplied by 10/ z s , and the r cc values by 



10 



-8.5-1/2 



We can see from the figures that the optical dis- 
tances initially increase rapidly with decreasing red- 
shift; this growth is slower for the absorption dis- 
tances at the frequency of the maximum of the Planck 
function. If the source's redshift is z 
z s /3 < z < z s , these distances are 



z s , then, for 



pi c z s - z 



T ~ (1 - e- CW ) ( -^-] 



h 



For the adopted model parameters, tt &0.13(z s — z), 
and t cc ~ 1.8 x 10~ 9 z<r 1/ ' 2 (;z s — z). The optical dis- 
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tances increase more slowly with decreasing z, gradu- 
ally reaching their asymptotic values tt ~ 0.14z s and 
t cc <~ 2.5 x lO -9 ^ 2 . The distances r cc reach their 
limiting values at smaller z, and these values are sub- 
stantially smaller (by a factor of 5.6 x lO 6 ^ 2 ) than 
those for the scattering distances. Since the coeffi- 
cient k cc is inversely proportional to the third power 
of the frequency this difference will be smaller at 
lower frequencies. 

4.7 The Grey Approximation 

Since the scattering coefficient substantially exceeds 
the absorption coefficient, the latter can be neglected 
relative to the source radiation. Equation (24) with 
the right-hand side (26) can then be integrated over 
frequency, as is done in calculations of model atmo- 
spheres [4]. 

Let us denote the frequency moments for the mean 
occupation number (s > is not necessarily integer, 
so that this is essentially the Mellin transform): 



n 



(«) 



OO 

{V, X, M) = J x s n(i], x, n, x)dx. 



After separating out the direct radiation, the function 
e* should be substituted by e s . 

4.8 CONCLUSION 

We have derived an equation describing the evolution 
of the radiation of a source whose intensity exceeds 
that of the surrounding equilibrium background ra- 
diation at the epoch of radiation-dominated plasma, 
between the epochs of annihilation and recombina- 
tion. A cosmological model for this period was 
adopted. We have separated the diffuse and direct 
radiation, and calculated the opacities of the mat- 
ter due to Thomson scattering and bremsstrahlung 
absorption. The latter mechanism can be neglected 
relative to the source radiation, and this fact was used 
to obtain an equation for the frequency moments. 

This is the first paper in the series of studies that 
will be concerned with elucidating how sources of ra- 
diation affect the evolution of the Universe and how 
their presence is reflected in the thermal background 
radiation, taking into account the fact that the accu- 
racy of observations is growing steadily. We hope to 
solve this problem in our forthcoming studies. 



Integration over the frequency with the weights x s 
yields 



Q n (,s) Q n (s) 



dr] 



dx 



srifc(x) djj, 



+(s + 1) 



a'(v) . 



= R(r,)4 s \ 



(s) 

where I c is the result of integrating the right-hand 
side. The integral of the term with the frequency 
derivative was calculated via integration by parts. 
The frequency moments of the collision integral are 



= - 



CUT 



n {s) {q,x,IJ) -n}?>(r),x) 



-^P2^)4 s \v,x) 



+ £ { * S \v,X,v), 



where the moments of the source power are 

oo 

el s) (»?, X, /") = J e* (V, X, Mi x)x s dx. 
o 
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